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Abstract 

Western lowland gorillas (Gorilla gorilla gorilla) are designated as critically 
endangered and wild populations are dramatically declining as a result of habi- 
tat destruction, fragmentation, diseases (e.g., Ebola) and the illegal bushmeat 
trade. As wild populations continue to decline, the genetic management of the 
North American captive western lowland gorilla population will be an impor- 
tant component of the long-term conservation of the species. We genotyped 26 
individuals from the North American captive gorilla collection at 11 autosomal 
microsatellite loci in order to compare levels of genetic diversity to wild popu- 
lations, investigate genetic signatures of a population bottleneck and identify 
the genetic structure of the captive-born population. Captive gorillas had signif- 
icantly higher levels of allelic diversity (f 7 = 4.49, P = 0.002) and heterozygosity 
(t 7 = 4.15, P = 0.004) than comparative wild populations, yet the population 
has lost significant allelic diversity while in captivity when compared to foun- 
ders (f 7 = 2.44, P = 0.04). Analyses suggested no genetic evidence for a 
population bottleneck of the captive population. Genetic structure results 
supported the management of North American captive gorillas as a single 
population. Our results highlight the utility of genetic management approaches 
for endangered nonhuman primate species. 



Introduction 

Wild populations of western lowland gorillas (Gorilla gor- 
illa gorilla) are facing dramatically declining numbers, as 
high as 80% over three generations, as a result of habitat 
destruction and fragmentation, the illegal bushmeat trade, 
and diseases such as Ebola (Walsh et al. 2008). Due to rap- 
idly increasing human encroachment into gorilla habitat, 
captive gorilla populations represent an important aid to 
long-term conservation as a representative species commu- 
nicating conservation issues in great ape range countries to 
the public. The captive collection of gorillas in North 
America began over 100 years ago with wild individuals 
imported from Africa; however, since coming under pro- 
tection of the Convention on International Trade in 
Endangered Species of Wild Fauna and Flora (CITES) in 
1975 there have been no new wild gorillas added to zoos 
(Nsubuga et al. 2010). Further, they were designated as 
critically endangered in 2007 (Walsh et al. 2008). 

In the wild, western gorillas (Gorilla gorilla diehli, 
G. g. gorilla) are separated from eastern gorillas (Gorilla 



beringei beringei, Gorilla beringei graueri) by the Congo 
River in central Africa. Eastern gorillas are listed as 
endangered at the species level and the mountain gorilla 
subspecies (G. b. beringei) is listed as critically endangered 
(Walsh et al. 2008). With the exception of the mountain 
gorilla populations in Bwindi and Karisoke, which have 
seen a population increase and stabilization in the last 
20 years due to "extreme" conservation efforts, other wild 
populations continue to decline (Guschanski et al. 2009; 
Robbins et al. 2011; Gray et al. 2007). Captive popula- 
tions, with no influx of wild individuals, may potentially 
face the same problems regarding loss of genetic diversity 
as small wild populations without appropriate genetic 
management (Ballou and Lacy 1995). The captive North 
American population of -370 individuals has the poten- 
tial to lose genetic diversity through random genetic drift 
if not properly managed as a single population, thus it is 
critical to have a breeding program that aims to maximize 
genetic diversity in order to avoid inbreeding depression. 
The gorilla Species Survival Plan (SSP) goals for 
maintaining genetic diversity are to maintain >90% of the 
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genetic diversity of the previous generation, over 
100 years, following Frankham et al. (2002). 

Breeding protocols of the North American Gorilla SSP 
include determining breeding pairs based on pedigree 
information (Nsubuga et al. 2010). However, accurate 
pedigree determination can be challenging because impor- 
tation records of wild caught gorillas can be inaccurate. 
For example, Nsubuga et al. (2010) discovered that a 
known breeding pair consisted of first order relatives 
from the founder population of wild-born captive gorillas. 
This is not an error of the zoos but reflects the difficulty 
of establishing kinship relationships in the absence of 
genetic data. The use of molecular markers, particularly 
microsatellite panels, has been shown to be a highly effec- 
tive genetic management tool in nonhuman primates 
(Deinard and Kidd 2000; Meier et al. 2000; Kanthaswamy 
et al. 2006; Perwitasari-Farajallah et al. 2010). Microsatel- 
lites are the marker of choice for genetic management for 
a number of reasons including high polymorphic infor- 
mation content and gene diversity. They are also effective 
at determining individual genetic identity and parentage 
exclusion (Kanthaswamy et al. 2006). These are particu- 
larly true for microsatellites with a tetranucleotide repeat 
motif, which are more reliably characterized than dinucle- 
otide repeats (Kanthaswamy et al. 2006). Genotyping 
individuals at polymorphic, neutral markers, such as 
microsatellites, and determining breeding pairs based on 
lowest mean kinship can help maximize genetic diversity 
in populations (Ballou and Lacy 1995). 

We investigated genetic diversity of the North American 
captive-born population of gorillas using a well-character- 
ized panel of microsatellite loci. We compared this new 
data set to previously collected data from wild gorillas as 
well as the wild-born founder population of captive North 
American gorillas to compare genetic diversity and investi- 
gate signatures of a genetic bottleneck. Because of high 
levels of genetic diversity and wide range of individuals in 
the wild-born founder population, we predicted that the 
genetic variation maintained in the captive-born popula- 
tion would be greater than that found in small wild 
populations of both western and eastern gorillas. Because 
of high levels of genetic diversity in the founder popula- 
tion, as a result of having multiple wild source 
populations (Nsubuga et al. 2010), we predicted a reduc- 
tion in genetic diversity in the captive-born population 
compared to the wild-born founder population. Lastly, 
based on the long generation time of gorillas and the 
relatively short time in captivity, we predicted captive- 
born gorillas, which are all western lowland gorillas, will 
cluster with wild western lowland gorillas (WWLG) to the 
exclusion of less related cross river gorillas and eastern 
gorillas, but that those clusters will reflect ancestral gene- 
flow between cross river and western lowland gorillas. 



Methods 

Samples 

DNA samples for captive western lowland gorillas (N = 26) 
were obtained from fibroblast cell lines located in the Inte- 
grated Primate Biomaterials and Information Resource col- 
lection at the Coriell Institute for Medical Research. 

Comparative datasets for this study came from five wild 
populations of gorillas and one founder population of the 
captive North American western lowland gorillas (Nsubuga 
et al. 2010). Included in the wild populations were Cross 
River gorillas (G. g. diehli) from the Cameroon-Nigeria 
border (Bergl 2006), and two populations of western 
lowland gorillas (G. g. gorilla) from Loango National Park, 
Gabon, and Mondika in the Central African Republic 
(Bradley et al. 2004; Bergl et al. 2008; Arandjelovic et al. 
2010, Fig. 1; respectively). The two western lowland gorilla 
populations were combined using the weighted mean of 
diversity measures and are referred to as wild western low- 
land gorillas. Eastern gorilla populations included two 
populations of mountain gorillas (G. b. beringei) from 
Bwindi Impenetrable Forest, Uganda and Volcanoes 
National Park in Virunga Mountains, Rwanda (Bradley 
et al. 2005; Nsubuga et al. 2008; Fig. 2). Based on sample 
sizes and sampling ranges we consider all four populations, 
and the founding population (WWLG, N= 131; Cross 
River, N = 71; Bwindi, N = 77; Virunga, N = 92; Foun- 
ders, N = 79) to be small, comparable populations. 

DNA amplification and genotyping 

Multilocus genotypes from 11 polymorphic autosomal 
microsatellite loci were produced for 26 individual goril- 
las. Microsatellite loci included Dls550, D2sl326, 
D5sl470, D4sl627, D5sl457, vWF, D16s2624, D8sll06, 
D10sl432, D2s2204, and D7s817 (Primers from Bradley 




Figure 1. Western lowland gorilla (Gorilla gorilla gorilla). Photo 
Credit: Nelson Ting. 
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Table 1. Allelic diversity and heterozygosity of captive-born gorillas. 



Cross River 



Mondika 



Virunga 



| Gorilla gorilla dtehli 
j Gorilla gorilla gorilla 
j Gorilla behngei beringei 
| Gorilla benngei grauen 



*Founder gorillas from Cameroon, Republic of Congo and Gabon 

Figure 2. Distribution map of populations included in the study. 
Adapted from Bergl et al. (2008) *Founder gorillas from Camerron, 
Republic of Congo and Gabon. 

et al. 2000). This panel of microsatellite loci was chosen 
based on the availability of comparable datasets in the lit- 
erature. All loci are tetranucleotide repeats with the 
exception of D5sl470, which is a tetranucleotide repeat 
with a 2 bp indel. 

PCR reactions were carried out with the following: 
10 /(L (5 U/^L) AmpliTaq Gold (Applied Biosystems, Fos- 
ter City, CA), 0.4 /iL (10 /imol/L) of forward and reverse 
primers, 7.7 fiL dH 2 0 and 1.5 fiL (10-30 ng//iL) DNA 
template for a final volume of 20 /.iL. Thermocycler proto- 
col included a 10 min initial denaturing at 95°C, followed 
by 40 cycles of 94°C denaturing for 30 sec, 55-60°C 
annealing for 30 sec (Table 1), 72°C extension for 30 sec, 
followed by a final extension of 72°C for 30 min. Forward 
primers were labeled at the 5' end with IR Dye 700 
fluorescent label. Amplicons were separated by gel electro- 
phoresis on a LiCor 4300 DNA Analyzer using 6.5% KB 
Plus gel matrix (LiCor Biosciences, Lincoln, NE). Alleles 
were scored using an IR Dye 700 internal size standard in 
SagaGT genotyping software (LiCor Biosciences). Geno- 
type data was checked for allelic dropout and null alleles 
using CERVUS 3.0.3 (Marshall et al. 1998; Kalinowski 
et al. 2007). 



Genetic diversity comparisons 

The captive-born gorilla dataset was examined for evi- 
dence of deviation from Hardy-Weinberg equilibrium 
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0.88 
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07.25 


7.03 


4.63 


0.86 


0.87 
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NA, number of alleles; AR, allelic richness; AE, number of effective 
alleles; H\, individual heterozygosity; H Q , observed heterozygosity; H E , 
Nei's expected heterozygosity. 



and linkage disequilibrium in GENEPOP 4.0.9 (Rousset 
2008). Genetic variation in the captive-born population 
was compared to the wild-born population of founder 
gorillas and four populations of wild gorillas using mea- 
sures of allelic diversity and heterozygosity. Of the 11 
micro satellites used, a subset of eight loci was used for 
inter-population comparisons based on the availability 
of comparable datasets in the literature. The number of 
alleles (NA), and number of effective alleles (AE; Kimura 
and Crow 1964), were calculated for the captive-born 
and Loango dataset (Arandjelovic et al. 2010) in GenAl- 
Ex 6.01 (Peakall and Smouse 2006). Number of effective 
alleles is a measure of the evenness of the allele 
frequency distribution averaged over all loci. This mea- 
sure is suited for comparing populations with differing 
numbers of alleles. Allelic richness (AR), a measure of 
alleles which controls for differences in sample size, was 
calculated for the captive-born and Loango populations 
in FSTAT 2.9.3.2 (Goudet 2001). Captive-born and wild 
populations were evaluated for differences in allelic 
diversity (NA, AR, AE) using t-tests. Levels of observed 
heterozygosity (H Q ), and expected heterozygosity (H e ), 
were calculated in GenAlEx 6.01 (Peakall and Smouse 
2006). Mean individual heterozygosity (H{), was calcu- 
lated similar to Nsubuga et al. (2008), as the mean 
number of heterozygous loci for each gorilla, divided by 
the total number of loci. Levels of heterozygosity in the 
captive-born population were tested for differences 
among the founding population and wild populations 
using f-tests following Bergl et al. 2008; Archie 1985; 
Nei 1987). Statistical significance was set at a = 0.05 
and corrected for multiple comparisons using the 
Holm's sequential Bonferroni adjustment (Holm 1979; 
Jaccard 1998). The Holm's sequential method is a more 
sophisticated correction that controls for inflation of the 
Type I error rate while also maintaining statistical power 
(Kromery and Dickson 1995). Adjusted levels of alpha 
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Table 2. Mean measures of allelic diversity and heterozygosity. 



Population 


NA 


AR 


AE 


Hi 


Ho 


H E 


Captive 


7.25 


7.03 


4.63 


0.86 


0.87 


0.76 


Founder 1 


8.50 


8.48 


5.18 


0.72 


0.73* 


0.80 


Wild western 


6.62 


5.82* 


3.88 


0.78 


0.75* 


0.68 


Lowland gorillas 2,3 














Cross River 4 


6.00 


4.91** 


3.30* 


0.68 


0.68** 


0.68 


Bwindi 5 


6.13 


5.06** 


3.29* 


0.71 


0.70** 


0.68 


Virungas 6 


5.13* 


4.22** 


2.69* 


0.71 


0.68** 


0.62* 



NA, number of alleles; AR, allelic richness; AE, number of effective 

alleles; H t , individual heterozygosity; H Q , observed heterozygosity; H E , 

Nei's expected heterozygosity. Significant values shown in bold. 

1 Nsubuga et al. (2010). 

2 Bradley et al. (2004), Bradley et al. (2007). 

3 Arandjelovic et al. (2010). 

4 Bergl et al. (2008). 

5 Nsubuga et al. (2008). 

5 Bradley et al. (2005). 

*P< 0.05. 

**P < 0.005. 

ranged from 0.003 to 0.05. The statistical significance of 
stated P values are relative to adjusted levels of alpha. 



Genetic structure and population 
assignment 

Genetic population structure in the captive-born popula- 
tion of gorillas was inferred using STRUCTURE 2.3.3 
(Pritchard et al. 2000). For structure analysis, a subset of 
five loci was used (D8sll06, D16s2624, vWF, D5sl457, 
Dls550) based on availability of genotype datasets for com- 
parative populations (Loango, Cross River, Virungas). This 
subset was used due to the Loango dataset containing geno- 
types from nested primers for five loci, which were 
removed. Due to differences is base-calling, the Loango 
dataset appeared to be called 2 bp below the other three 
populations for the remaining five loci. When the Loango 
dataset was corrected for this the allele frequency distribu- 
tion matched for all four populations and the datasets were 
compatible. The program STRUCTURE uses a Bayesian 
model-based clustering method to infer population genetic 



structure under the assumption of K clusters, where K is 
the number of individual clusters or populations. In order 
to determine the optimal number of K for this dataset, 10 
independent replicates were run for values of K = 1-5. For 
each replicate within each K value, tests were run with 
100,000 burn-in steps followed by 1,000,000 Markov Chain 
Monte Carlo replications. All tests were run under an 
admixture model with correlated allele frequencies. The 
log-likelihood [In P(D)] was averaged over independent 
runs for each value of K. Runs of high K values can poten- 
tially increase the posterior probability as well as variance 
between independent runs, leading to an overestimation of 
K (Rosenberg et al. 2001; Nsubuga et al. 2008). Following 
the method of Evanno et al. (2005), we used the ad hoc test 
statistic AK, which is the second order rate of change in In 
P(D) across consecutive values of K. The use of AKto iden- 
tify breakpoints in the dataset results in the true value of K 
being that with the greatest AK. 

Population assignment tests were carried out with the 
four datasets used for inferring genetic population 
structure (captive-born, Loango, Cross River, Virungas) in 
order to assess whether the captive-born population would 
cluster with the Loango population, which is the closest 
population to founders as well as being the same sub-spe- 
cies, to the exclusion of the Cross River and Virungas. 
Assignment tests were carried out in GenAlEx 6.01 (Peak- 
all and Smouse 2006). Pairwise population assignment 
tests were also conducted among captive-born population 
and Cross River, Loango and Virunga populations. Assign- 
ment tests were frequency-based following Paetkau et al. 
(1995). Assignments tests were based on the log-likelihood 
value of genotype frequency over all loci for each popula- 
tion. Individuals are then assigned to the population with 
highest log-likelihood. 

Demographic history 

Signatures of a genetic bottleneck were tested in the 
captive-born population using the BOTTLENECK 
program (Piry et al. 1999). Because the mutation model 
for these microsatellite loci is not known, we used a 



Figure 3. Structure results including captive- 
born and wild populations of western lowland 
gorillas, Cross-River gorillas, and mountain 
gorillas. 
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two-phase mutation model (TPM), which combines the 
infinite alleles model (IAM) and step-wise mutation 
model (SMM). The TPM accounts for the unlikelihood 
that microsatellite loci will precisely follow either a strict 
SMM or IAM (Di Rienzo et al. 1994; Piry et al. 1999). 
Following the method of Weckworth et al. (2005) we per- 
formed runs with step-wise changes in the contribution 
of SMM to the TPM of 70%, 75%, 80%, 85% and 90%. 
Additionally, we ran the test under a strict IAM and 
SMM separately. Three tests were used to assess signifi- 
cance in the difference between H E and H EQ , where H E is 
the expected heterozygosity, assuming mutation- drift 
equilibrium, and H EQ is a coalescent-based estimate of 
heterozygosity based on the observed NA (Piry et al. 
1999). Tests used for significant differences in H EQ and 
H E were Wilcoxon sign-rank tests, sign test and standard- 
ized differences test. 

Results 

Genetic diversity in captive-born gorillas 

No evidence of null alleles or allelic dropout was observed 
in the captive-born dataset. In both the captive-born goril- 
las (D4sl627) and Loango (vWF), a single locus was 
observed to deviate from Hardy-Weinberg equilibrium. 
Deviations that have been previously described in Cross 
River (D5sl470 and D8sll06), Bwindi and Virunga 
populations (Dls550 and D4sl627, respectively) were 
attributed to the inclusion of related individuals in the 
sample (Bergl et al. 2008). In those datasets it was shown 
that when closely related individuals were removed, those 
loci no longer deviated from equilibrium (Lukas et al. 
2004; Bradley et al. 2005). Because the captive-born data- 
set includes closely related individuals, including three 
pairs of full siblings and two pairs of half siblings, we 
followed Bergl et al. (2008) in treating all loci as though 
they were in equilibrium and independent. 

The captive-born population showed relatively high 
levels of both allelic diversity and heterozygosity measures 
(Table 1). The captive-born gorilla population had higher 
levels of allelic diversity than that found in wild popula- 
tions (Table 2). The NA was significantly higher in 
captive-born gorillas than in Virungas. For measures of 
AR, the captive-born gorillas were significantly higher 
than Virungas, Bwindi, Cross River and WWLG popula- 
tions. Cross River and Virungas populations had 
significantly lower AE than captive-born. Captive-born 
gorillas had lower levels of AR than the wild-born foun- 
ders. Captive-born AR was significantly lower than the 
wild born founders but NA was not significantly differ- 
ent. Effective alleles in captive born gorillas were also not 
significantly lower than founders. 



Population assignment 
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Figure 4. Frequency based population assignment of four 
populations including captive-born, Cross River, Loango and Virunga 
(Karisoke) gorillas. 



Population assignment for captive vs. loango 
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Figure 5. Pairwise frequency based population assignment plot 
including captive-born and Loango gorillas. 



Captive-born gorillas had higher levels of two hetero- 
zygosity measures than the wild-born founders. The H Q 
was higher in the captive-born gorillas than in the foun- 
ders. Captive-born gorillas also had significantly higher H Q 
than wild populations. Similarly, H e was significantly 
higher in the captive-born gorillas than in Virungas. 
Captive-born gorillas did not significantly differ from the 
founders in H e . 

Population structure and assignment 

Comparisons of the second order rate of change, AK 
(Fig. 3) of the In P(D) from STRUCTURE found a major 
breakpoint in the data with the highest likelihood of clus- 
ters at K = 3 (Pritchard et al. 2000; Evanno et al. 2005). 
The captive-born population clustered with the Loango 
population with the majority of individuals having a pro- 
portional group membership value (Q) > 80% (Fig. 3) at 
K = 3. Because captive gorillas and Loango gorillas are 
the same subspecies, we expected them to cluster 
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together, which resulted in a K = 3, when combined with 
Cross River and Virunga populations. 

Frequency based assignment tests agreed with the 
Structure results and clustered captive-born and Loango 
populations together to the exclusion of Cross River and 
Virungas (Fig. 4) when all four populations were 
included. Pairwise population assignments showed cap- 
tive-born and Loango gorillas formed less discrete clusters 
than captive-born and other populations (Figs. 5-7). As 
the source for the captive populations are thought to 
come from Cameroon, the Congo and Gabon, we would 
expect that captive-born and Loango gorillas would clus- 
ter more closely than other populations to captive-born, 
which they do in the pairwise assignment, suggesting the 
methodology and markers used are accurate in assigning 
individuals to these populations. 

Genetic bottleneck signatures 

The results of the BOTTLENECK analysis differed based 
on model choice but did not suggest the occurance of a 
bottleneck. The BOTTLENECK results using the TPM 
model were not significant according to any of the three 
tests, under any of the step-wise runs. Similarly, under a 
strict SMM, all three tests were not significant. In con- 
trast, under a strict IAM, the Wilcoxon sign-rank test was 
significant for a one-tailed test of H EQ excess (P = 0.002). 
The standardized differences test was also significant 
(P = 0.005). These results were not surprising as only a 
severe bottleneck would be detectable after only two gen- 
erations. Further, if there were evidence of a bottleneck it 
would more likely reflect the demographic history of the 
source populations. 

Discussion 

Comparative genetic diversity of the 
captive-born population 

Concerning the genetic diversity of the North American 
captive population of western lowland gorillas, they show 
high levels of both allelic diversity and heterozygosity. 
The captive-born population had a higher mean NA, AR 
and effective alleles than the wild populations. The cap- 
tive-born population was significantly higher than all wild 
populations in mean AR, which is the most informative 
measure between these datasets. This meets the expecta- 
tion that the captive-born gorillas would have high levels 
of allelic diversity when considering that the founder 
gorillas also have higher measures of allelic diversity than 
all wild populations. For all three measures, NA, AR and 
AE, the founder population had higher values than the 
captive-born, suggesting that while still high, the captive 



Population assignment for captive vs. cross river 




Figure 6. Pairwise frequency based population assignment plot 
including captive-born and Cross River gorillas. 

gorilla population has lost some allelic diversity present 
in the founder population. H Q was significantly higher in 
captive-born than founders, despite there being no signifi- 
cant difference in H e . This observation of decreased allelic 
diversity but increased H a from the founder to captive- 
born population is interesting. We interpret this observa- 
tion to be an artifact of changes to breeding protocols. 
Until recently breeding pairs were determined based on 
flawed pedigrees and closely related individuals were 
known to have bred (Nsubuga et al. 2010). This could 
have resulted in the loss of rare alleles, which may 
account for the decrease in allelic diversity from the foun- 
der to captive-born population. Recently, breeding based 
on non-random negative assortative mating according to 
least mean kinship between pairs would result in an 
increase in H Q from the founders, while H e (which is 
based on NA) and NA have not had enough time to 
rebound through mutation. For measures of both allelic 
diversity and heterozygosity, the captive population is 
higher than all wild populations, despite the having lost 
allelic diversity since being in captivity. 

Genetic structure and population 
assignment 

Based on the results of the STRUCTURE analysis, cap- 
tive-born gorillas form two clusters, which is consistent 
with the clustering of the wild-born founders (Nsubuga 
et al. 2010). Structure results also found that when 
grouped with the Cross River, Loango and Virunga popu- 
lations, the captive-born gorillas clustered with the 
Loango gorillas for a highest likelihood of three clusters. 
Assignment tests agreed with structure results that the 
captive-born and Loango gorillas were clustered together 
to the exclusion of Cross River and Virunga gorillas, and 
pairwise assignment test showed Loango and captive-born 
gorilla clustering more closely than captive-born and 
Cross-River or Virunga. 
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Population assignment for captive vs. virunga 
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Figure 7. Pairwise frequency based population assignment plot 
including captive-born and Virunga gorillas. 



Genetic bottleneck signatures 

All tests under both a TPM and strict SMM were not sig- 
nificant. Results from both the Wilcoxon sign-rank test 
and standardized differences test under a strict IAM were 
significant, yet the IAM is the least conservative mutation 
model. The SMM, and to a greater degree, the TPM are 
considered more informative models in addressing ques- 
tions of demographic history. These results suggest that 
there is no genetic evidence for a population bottleneck 
in this captive collection. While the population has lost 
allelic diversity, enough diversity has been maintained 
in the population that they are still significantly higher 
than comparable wild populations for several diversity 
measures. These results also suggest that despite flaws in 
pedigree information on founder gorillas, genetic manage- 
ment of this population has been successful in maintain- 
ing high levels of genetic diversity. 

Implications for conservation 

There is currently no scientific justification for the use of 
the North American captive gorilla population as a source 
for reintroduction of individuals to the wild (Beck et al. 
2007). As such, there is currently only one captive breed- 
ing program outside of a range country that re-introduces 
captive-born gorillas to the wild. The goal of the gorilla 
SSP is to manage the captive population of North Ameri- 
can gorillas as single unit, and our results confirm those of 
Nsubuga et al. (2010) that the captive population should 
be managed as a single population. Measures of genetic 
diversity in the captive gorilla population are high com- 
pared to other captive mammal populations. The data 
presented here may be useful in the continued monitoring 
of genetic diversity in the captive population for the 
long-term maintenance of zoo collections. The data fur- 
ther suggests that the set of eight microsatellite loci used 



here represent an ideal panel for genetic management use; 
they are highly polymorphic, relatively neutral when 
accounting for closely related individuals in the sample 
and are all tetranucleotide repeat motif. In addition, as 
presented here, there exists a number of comparative wild 
datasets using this panel of loci, which can be used to esti- 
mate changes in genetic diversity over time in the captive 
population relative to wild populations. 

Data Accessibility 

• File with genotype dataset. 

• Readme.txt file for above referenced data file. 
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